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We present a field theory for the statistics of charge and current fluctuations in diffusive 
systems. The cumulant generating function is given by the saddle-point solution for the 
action of this field theory. The action depends on two parameters only: the local diffusion 
and noise coefficients, which naturally leads to the universality of the transport statistics for a 
wide class of multi-dimensional diffusive models. Our theory can be applied to semi-classical 
mesoscopic systems, as well as beyond mesoscopic physics. 



1 Introduction 

In this short paper, we present the essential part of our theory^ for the statistics of current and 
density fluctuations in non-equilibrium diffusive systems. The theory is based on the stochastic 
path integral approach to the statistics of fluctuations in networks 1 ^ and represents an alternative 
to quantum method for the evaluation of the full counting statistics (FCS) of the transmitted 
charge, as well as to classical cascade correction method^for the evaluation of current cumulants. 

The building blocks of our theory are a separation of the time scales (sources of noise 
are Markovian) and a saddle-point approximation of the resulting functional integral. The 
separation of variables into fast currents and slow varying generalized charges suggests the 
network description of a system introduced in Sec. El In Sec. [3 we consider a large number of 
nodes and derive the continuum limit for the stochastic path integral, which turns out to be 
a stochastic field theory with an action being a function of only two quantities, the diffusion 
D and noise F functionals of the charge density p. This naturally leads to the universality of 
the FCS not only for the mesoscopic diffusive conductor [see Eq. (|22|)]. but for a whole class of 
multi-dimensional diffusive models (see also Ref. [1] for details). 

We would like to stress that due to its essentially classical construction, our stochastic 
field theory is not limited to mesoscopic physics^ but can also be applied to reaction-diffusion 
systems^ symmetric exclusion processes^ and many other areas of classical stochastic processes^ 

2 Stochastic Path Integral for a Network 

Consider a network with the state of each node a described by one charge Q a (Q is the charge 
vector describing the charge state of the network). The node's state may be changed by trans- 
port: flow of charges between nodes takes place via the connectors carrying currents I a p from 
node a to node j3. The rate of change of these charges Q a is given by 



where P a p is the probability distribution of the fast current between nodes a and (3, and H a p 
is the current cumulant generating function. The fact that P a p also depends on the charges Q 
is one source of the difficulty of the problem. 

Assuming that the generators H a p of the fluctuating currents I a p are known, we seek an 
evolution of the probability distribution T(Q,t) of the set of charges Q for a given initial con- 
dition r(Q,0). In other words, one has to find the conditional probability (which we refer to as 
the evolution operator) C/(Q,Q',t) such that 
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We assume that there is a separation of time scales, tq <C to, between the correlation time of 
current fluctuations, ro, and the slow relaxation time of charges in the nodes, re*. In Ref. [2] we 
have used the separation of time scales to derive a stochastic path integral representation for 
the evolution operator, 

U{Q f ,Qi,t) = |pQPAexp{5(Q,A)}, (3) 
5(Q,A) = f dt'[-iA -^ + (1/2)^2 H a p(Q,X a -X p )]. (4) 
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The variables X a are auxiliary variables for every node that impose charge conservation in the 
network. 

3 Continuum Limit 

From the stochastic network, it is straightforward to go to spatially continuous systems as the 
spacing between the nodes is taken to zero. Consider a series of identical, equidistant nodes 
separated by a distance Az. Th is nodal chain could represent a chain of chaotic cavities, Fig. 
El in a mesoscopic context P^-^ The sum over a and becomes a sum over each node in space 
connected to its neighbors. The action for this arrangement is 

S= [ d£ V]{— AqQq: + H(Q a , Q a -i; X a — A a _i)} , (5) 

where for simplicity we have chosen real counting variables, iX a — » X a . The only constraint 
made on H is that probability is conserved, H(X a — A Q _i) = for X a = X a -i. We now derive 
a lattice field theory by formally expanding H in X a — X a -i and Q a — Q a -\. Only differences 
of the counting variables will appear in the series expansion, while we must keep the full Q 
dependence of the Hamiltonian. If there are N S> 1 nodes in the lattice, for fixed boundary 
conditions the difference between adjacent variables, X a — X a -i and Q Q — Q a -i will be of order 
1/iV, and therefore provides a good expansion parameter. The expansion of the Hamiltonian 
© to second order in the difference variables gives 

H = -qj— (A a — X a -l) + - ~Q)2~(^a — A a -l) 2 + QQ-Q^iQot — Qa-l)(X a — A a _l) , (6) 

where the expansion coefficients are evaluated at A a = A a _i and Q a = Q a -i and are functions 
of Qa-i- Terms involving only differences of Q a — Q a -i are zero because H(X a — A a _i) = for 
X a = X a -i. All terms in Eq. (JHJ) need explanation. First, the expression dH/dX a is the local 
current at zero bias (because the charges in adjacent nodes are equal) which will usually be zero. 
There may be exceptional circumstances where this term should be kept, but we do not consider 
them here. The term <9 2 H/dQ a dX a = —G(Q a -i) is th e lin ear response of the current to a 
charge difference. Hence, G is the generalized conductance^^ of the connector between nodes a 
and a — 1. d 2 H/dX^ = C(Q a —i) is the current noise through the same connector because H is 
the generator of current cumulants. 

We are now in a position to take the continuum limit by replacing the node index a with a 
coordinate z, introducing the fields Q(z), X(z), and making the expansions 

X a - X a ^^X'Az + (l/2)A"(Az) 2 + 0(Azf, (7) 
Q a - Qa-i^Q'Az + (l/2)Q"(Az) 2 + 0{Azf. (8) 

The action may now be written in terms of intensive fields by scaling away Az, 

H^h(p,X)Az, Q a ^p(z)Az, G a (Az) 2 — > D(p), C a Az -> F(p) , (9) 




Figure 1: A one dimensional lattice of nodes connected on both ends to reservoirs. This situation could represent 
a series of mesoscopic chaotic cavities connected by point contacts. 



and taking the limit — > J dzh(p, A). One may check that expanding the Hamiltonian 
to higher than second order in Az will result in terms suppressed by powers of Az/L and 
consequently vanish as Az — » 0. 

These considerations leave the one dimensional action as 



S = - f dt' [ L dz 
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Xp + Dp'X' - l -F{\'f 



(10) 



Here D is the local diffusion constant and F is the local noise density which are discussed in detail 
below. It is very important that these two functionals D,F are all that is needed to calculate 
current statistics. Classical field equations may be obtained by taking functional derivatives of 
the action with respect to the charge and counting fields: 

A = -~ (A') 2 - DX", p = [-FX' + Dp'}' . (11) 

We have to solve these coupled differential equations subject to the boundary conditions p(t, 0) = 
PL(t), p(t,L) = pn(t), X(t,0) = Ax(i), and X(t,L) = Xn(t). Functions A^(i) and A_r(£) are the 
counting variables of the absorbed charges at the left and right end of the system. Once Eqs. (|llj) 
are solved, the solutions p(z,t) and X(z,t) should be substituted back into the action lfTU|) and 
integrated over time and space. The resulting function, S sp [pL(t), pn(t), Xi(t), XR(t),t, L] is 
the generating function for time-dependent cumulants of the current distribution. Often, the 
relevant experimental quantities are the stationary cumulants. These are given by neglecting the 
time dependence, finding static solutions, p = X = 0, and imposing static boundary conditions. 
We can also introduce sources J dtdz x( z > t) p(z , t) and calculate density correlation functions. 

To justify the saddle-point approximation, it is useful to define dimensionless variables. The 
boundary conditions pl, and pr provide the charge density scale po in the problem, so we define 
p{z) = pof(z), where / ~ 1 is an occupation. We furthermore rescale z — > Lz, and t — > Tpi, 
where td = L 2 /D is the diffusion time, thus obtaining 



S = -Lp f V C dz' 
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A/ + /A'--^(A') 2 



(12) 



We assume that the combination F/Dpo is of order 1. From Eq. (|12|). the dimensionless large 
parameter is 7 = pqL ^> 1, i.e. the number of transporting charge carriers. The saddle-point 
contribution is of order jt/ro, while the fluctuation contribution is of order t/ro- 
Repeating these steps in multiple dimensions yields the action 

S = -f t dt'( dr [Xp + VXDVp- (1/2)VAFVA] , (13) 
Jo Jn 



where F and D are general matrix functions of the density p and coordinate r which should be 
interpreted as noise and diffusion matrices. 



As in any field theory, symmetries of the action play an important role because they lead 
to conserved quantities. We first note that the Hamiltonian h(p, Vp, VA) is a functional of VA 
alone with no A dependence. This symmetry is analogous to gauge invariance, and leads to the 
equation of motion 

P + V-j = 0, j = -DVp + FVA , (14) 

which can be interpreted as conservation of the conditional current density j. The next sym- 
metry is related to the invariance under a shift in the space and time coordinates {<5r, St}. 
This symmetry leads to equations analogous to the conservation of the local energy/momentum 
tensorM^For the stationary limit (where p and A vanish) and for symmetric diffusion and noise 
tensors, the conservation law is relativly simple and is given by 

X] ^mTmn = , T mn = j m (V„A) - (V n p) (DVX) m - h 5 mn . (15) 

m 

For the special case of a one dimensional geometry, the Hamiltonian itself is the conserved 
quantity (see Sec. 0J). 



4 FCS of Diffusive Systems 

We first consider the general ID field theory with the action 1)10(1 . and then demonstrating 
our solution for the FCS of the mesoscopic diffusive wire specifically. In the stationary limit, 
p = A = 0, the action can be written as 



S = t J dz 
o 

The stationary saddle-point equations 

(FX' -Dp' 



-Dp'X' + l -F{X'f 



(16) 



0, 2DX" + 5 -f{X'f 
5p 



0, 



(17) 



can be partially integrated leading to the following two equations: 

Dp' = ±Vl 2 - 2HF, (18) 

X' = 2H/(1 - Dp'). (19) 

The two integration constants X = —Dp' + FA' and Ti = —Dp'X' + (F/2)(X') 2 are the conserved 
(conditional) current and the Hamiltonian density, respectively. These conservation laws fol- 
low from the symmetries of our ID field theory [see Eqs. (|14|) and (|15|) and the surrounding 
discussion] . Thus we obtain the following result for the action (TTH|) , 



S = tLH. (20) 

The Eqs. (|18H20(I represent the formal solution of the FCS problem for ID diffusion models 
with D(p) and F(p) being arbitrary functions of p. The following procedure has to be done 
in order to obtain the cumulant generating function S(x) of the transmitted charge: (i) The 
differential equation (|T8|) has to be solved for p(z) with the boundary conditions p(z)\ z= o = p^ 
and p(z)\ z= L = Pr- The constant X should be expressed through the constants pi, pr, and 
Ti. (ii) Next, p(z) is substituted into Eq. (|19|) which is integrated to obtain X(z) with the 
boundary conditions Xl = and Xr = \- (hi) Finally, using the solution for X(z) the constant 
Ti is expressed in terms of pi, Pr, X, an d substituted into the action (f2"U|) , We note that by 
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Figure 2: The logarithm of the distribution of the current through a mesoscopic diffusive wire as a function of 
the ratio I/Io of the current to its average value 7o. The distribution is strongly asymmetric, with the Gaussian 
tale at I ^> Iq. Inset: The electron occupation / inside the wire as a function of the rescaled coordinate z, under 
the condition that the average current I — Io, no current 1 = 0, and large current / = 5/o has been measured. 

expressing 7i and x m terms of X, we may also formally obtain the logarithm of the current 
distribution, 

lnP(I)=S(2)-tl X (l), !->/, (21) 

as a result of the stationary phase approximation for the integral P(J) = / d X exp[S(x) — tlx] 
and because dTL/dx = T/L. 

As an example of the ID field theory, we consider the FCS of the electron charge transmitted 
through the mesoscopic diffusive wire. When the potential difference A/j, = (jll — fJ-R > is 
applied to the wire, the electrons flow from the left lead to the right lead with the average 
current Iq = e _1 GA/i, where G is the conductance of the wire. The elastic electron scattering 
causes non-equilibrium fluctuations of the current. At zero temperature, and for noninteracting 
electrons (the cold electron regime), the FCS of the transmitted charge has been studied in 
Refs. [14] and [4] using quantum-mechanical methods with the following result for the generating 
function of cumulants of the dimensionless charge Q/e: 



S(x) = (tlo/e) arcsinh 2 



\J exp(x) 



(22) 



Here we will rederive this result using our classical method. 

On the classical level, the electrons in the diffusive wire are described by the distribution 
function f(z). Under transport conditions (and at zero temperature), this distribution f(z) 
varies from //, = 1 in the left l ead to fn = in the right lead. Taking the continuum limit for 
the series of mesoscopic cavitiesp^ we arrive at the action Ijlfijl in the form 

1/2 

S = (tl /e) J dz[-f'\' + f(l-f)(\')% (23) 

-1/2 

where we have rescaled the coordinate z, p(z) has been replaced with the distribution f(z), and 
where D = 1, and F = 2/(1 — /) up to the overall constant I$/e. This form of F, originating from 
the Pauli blocking factors, is quite general for fermionic systems. Applying now the procedure 
described in the beginning of this section, we solve the saddle-point equations and find the fields 



/ and A, 



1 



sinh(2az) 



sinha 

X(z, x) = 2 arctanh [tanh(a/2) tanh(az)] 



a = arcsinh 



\/exp(x) 



1 



(24) 
(25) 
(26) 



where TC = a 2 , so that according to the Eq. (|20j) we immediately obtain the result (I22j) . 

The logarithm of the current distribution ln[P(/)] can be now found from the equation (J3 
We obtain the following result: 



a 



(27) 



ln[P(/)] = -(t/ /e)[2acothaln(cosha) 

where a has to be expressed in terms of 1 = I/Iq by solving the equation 

acotha = I/I . (28) 

The distribution P(I) is strongly asymmetric around the average current I = Iq (see Fig. |2J). 

In Ref. [1] we have proven the universality of the FCS of the transmitted charge for a two- 
terminal multi-dimensional generalized wire with the noise tensor F(p)T, being an arbitrary 
function of the charge density p, and with the constant diffusion tensor DT. The universality 
means that the FCS depends neither on the shape of the conductor, nor on its dimensionality. 
The FCS of a mesoscopic wire given by Eq. (|22|) is a particular example of universal FCS. In 
the more general case, when D is a function of p, the FCS depends on the geometry through 
only one parameter, the geometrical conductance. 
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